Attosecond streaking of core lines of copper dihalides 
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In the attosecond (as) streaking of Cu 3s core-level photoemission of copper dihalides, we predict 
theoretically that the satellite (3d^) is emitted later than the main line (3d^°L-^ L: ligand). The 
emission time delay is originated from the electron correlation between the core level and 3d shell, 
which leads to the difference in core-hole screening between satellite and main hnes. Further, we 
find that the time delay corresponds to a quantification of the extrinsic loss of photoemission. 
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During the past decade, advances in attosecond (as) (1 
as = 10^^^ s) laser technology have allowed time-resolved 
tracking at a time scale of the order of 0(10) - 0(100) 
as of the photoinduced electron dynamics occurring in 
atoms or molecules ^ perspective has been re- 
cently opened up by attosecond streaking metrology with 
the photoemission from the condensed phaseQ. Photoe- 
mission spectroscopy is one of the most developed tools 
for investigating the electronic structure of condensed 
matter 0. The combination of attosecond metrology and 
photoemission spectroscopy enables time-domain insight 
into the electron photocreation and transport stages of 
the photoemission process. 

In attosecond streaking photoemission, an electron 
emitted by an attosecond extreme ultraviolet (XUV) 
pulse is streaked by the femtosecond infrared (IR) pulse 
of a variable relative delay. The first proof of the prin- 
ciple experiment has revealed a time delay of about 100 
as (100 ± 70 as) between photoelectrons from the local- 
ized 4/ core level and the conduction band of a tung- 
sten surface^. A large part of the observed time delay 
is found to be from the transport delay due to a dif- 
ference in group velocities of electrons emitted from the 
core level and conduction band 5, 7, 8]. In a subsequent 
attosecond photoemission experiment for the atomic tar- 
get, a different kind of time delay was observed. The 
slower 2s electrons are found to be emitted earlier than 
the faster 2p electrons by 21 ± 5 in neon atoms [o^. This 
time delay is not explained by the transport, but is re- 
lated to the emission timing [1, The energy deriva- 
tive of the quantum phase of the dipole matrix element 
determines this emission time delay, called the Wigner- 
Smith delay [H [12. The Wigner-Smith delay observed in 
neon atoms could provide a new insight into intra-atomic 
electron correlation, because it eventually determines the 
quantum phase. This naturally tempts a challenge to- 
ward a new time-domain understanding of electron cor- 
relation in a true solid with strongly correlated electrons. 

Copper dihalides are insulating solids of id transition 
metal compounds, the simplest and best understood sys- 
tems where the charge transfer scenario is applied 13] . 
Figure [Ha) shows that photoexcitation of the Cu core 
level results in two final states. The Sd^ final state leaves 
the valence configuration largely unaltered and leads to 
a satellite (excited one), whereas the 3d^°L~^ final state 
(L: ligand shell) transfers an electron from a ligand to 




FIG. 1: (a) Copper dihahdes before and after photoemission. 
For instance, d implies a Cu 3ci state, c a Cu (3s) core level, 
and L a ligand, i.e., L — (3p)^ of CI for CuCb. U is the 
electron correlation between c and d. (b) Crystal structure of 
CuCl2: a CuCle cluster is designated by a blue dashed circle, 
(c) Vs,(r) for CuCla. 



a metal ion and leads to the main line (ground one) |l4| . 
Within this simple model, the extrinsic loss of photo- 
electron, that is, loss effects due to an interaction of 
photoelectron with the remaining solid during its travel 
to the surface, was studied in copper dihalides by Lee 
et a/. [15]. Their work could provide the key inspiration 
to the problem of attosecond streaking photoemission of 
strongly correlated electron systems. 

In this Letter, we study the attosecond streaking of 
3s core lines of copper dihalides by explicitly consider- 
ing the photoelectron scattering due to the core hole left 
after photoemission. We find that electron correlation 
between the core level and 3d shell leads to a difference 
in core-hole screening between the satellite and the main 
lines, which enables the theoretical prediction that the 
satellite is emitted later than the main line, say by 90 as 
(for CuBra), 73 as (for CuClj), and 40 as (for CuFa) at 
'^xuv — 141 eV, where wxuv is the XUV photon energy. 
The emission time delay strongly depends on the XUV 
photon energy and disappears at just ~ 20 eV over the 
threshold (i.e., the binding energy of the satellite). We 
further find that the observed time delay is a quantifica- 
tion of the extrinsic loss of photoemission. 

Photoemission analysis of strongly correlated electron 
systems including d electrons could be often simplified 
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from their localized nature. The simplest, but quite suc- 
cessful model based on a molecular orbital approach, rep- 
resents the d part of the system by one orbital, and the 
ligand by another jT3|. We start from a Hamiltonian Ho 
of three-state (two-state -I- one core level) model describ- 
ing copper dihalides, being consistently with FiglD^a), 

Ho = ednd+eLnL+ecnc + Uncnd+t{c]jCL + cl^Cd), (1) 

where the first two terms depict energy levels of the Cu 3d 
shell {ud = cj^Cd) and a halide ligand L {til = c]^cl) and 
the last term depicts the hybridization between them. 
The remaining terms involve the energy level of Cu 3s 
core level {nc = cj^cc) and the electron correlation U 
between the core level and the 3d shell. From Eq.(l), 
the initial ground state lipo) is readily obtained by in- 
serting nc = 1 as \ipo) — — single) Id) -I- COSMIC) |L) with 
tan 26* = 2t/{ed + U ~ El)- The ground state energy Eq 
is Eo=ec + ^{ed + C/ + sl) - ^^{ea + U ~ el)^ + 4t^. 
The final states of the target (after photoemission) are 
also immediately obtained by inserting nc = as 
hAm) — cosip\d) — sin ip\L) and |^s) — sin ip\d) + cos (p\L) 
with tan2(/j = 2t/(eL — £d), corresponding to the main 
{?>d^^L~^) and satellite state (3d^), respectively, as illus- 
trated in Figllja). The energy eigenvalues and Eg are 
iJ™ = \{ed + where 5E is the separation be- 

tween two states given by 5E — yj {sd — elY + Then 
photocurrent is J(w) — \wra\^5{ijj + Eo—E^-a) + \ws\^5{ijj + 
Eq — Eg), where Wm = — sin((/3 -I- 6) and Ws — cos{(p + 9) 
from u;,„(s) — (V'm(s) |cc|V'o)- Ji^) can be directly com- 
pared with the left panel of Figl5]so that the material pa- 
rameters would be taken from the relative peak strengths 
and peak positions in the experiment [l^ . This analysis 
constructs the sudden approximation. 

Hi comprises the excitation of photoelectron by the 
XUV pulse and its streaking by the IR pulse. 

Hi = ^ [ek - k • AiR(r -I- riR_xuv)] "k 

k 

+ rn]^{cl.cc + c|gCk)Axuv(T), (2) 

k 

where rik = cj^Ck is the occupation number of 
photoelectron with kinetic energy of £k — k^/2. 
AiR(r) is the streaking IR pulse given by Air(t) = 
0.2exp[-3.24 x lO^^r^] cos(0.0608T)eiR and Axuv(t) 
is the photoexciting XUV pulse given by Axuv(''') = 
A§;uvexp[-0.036T2]cos(wxuvT)exuv- XUV and IR 
pulses are assumed to be linearly polarized along the 
direction of photoelectron detection, tir and fxuv are 
the half width at half maximum (HWHMs) of the IR 
and XUV pulses, i.e., tir = 146.26 a.u. = 3.5 fs and 
^xuv — 4.39 a.u. = 105 as (where a.u. means atomic 
unit). TiR-xuv is the relative delay of the IR pulse from 
the XUV pulse, is (k|A|C) of the dipole operator A 
(= r • exuv) with |k) as the plane wave for the photo- 
electron, which is equivalent to the so-called strong field 
approximation. Noting that the k-dependence of Wk is 
usually weak, we take just as a constant. 
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FIG. 2: Left: Photoemission spectra of Cu 3s core lines of 
copper dihalides [T3]. Right: Centers of energy (COEs) of 
streaked photoemission spectra of copper dihaUdes with re- 
spect to a relative delay tir-xuv of the IR pulse from the 
XUV pulse (tJxuv = 141 eV). Main lines (red lines) are prior 
to the satellites (blue lines). 



Beyond the sudden approximation, we should consider 
the photoelectron scattering V by the core hole asfisj 
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where U^^k'(= / dr'0^(r)I/('') (r)V'k' (r)) is the matrix el- 
ement of the Coulomb potential from the orbital us- 
ing the plane wave basis. The second line is obtained 
by assuming that the Cu 3d level is so local as to 
have U("^)(r) « V'^^'^^r) and noting Ud + ul = 1- We 
now have V^^, = J drip^{r)Vsc{r)iJk>{r) and Vsc{r) = 
V^^^r) - y('*)(r). V^'^\r) for the Cu 3d orbital, which 
can be calculated by generating the charge densities in 
terms of Slater's or bitalsfTvl, while V^^\r) would be ap- 
proximated by a spherical charged shell made of six lig- 
and orbitals of a CuClg cluster (for CUCI2) with the ra- 
dius i?o as in Figjljb). Rq is an average distance be- 
tween Cu and CI. In the lowest approximation, there 
would be no scattering outside the cluster before reach- 
ing the surface unlike the metallic casefl^. This gives 



r), where &{Ro — r) 
Vscir) for CuCla is 



VUr)='^[-V^''Hr) + ^\eiRo 
is the Heaviside step function [Tf 
displayed in FigUJc). 

Under the total Hamiltonian of Ho + Hi+V, by solving 
the time-dependent Schrodinger equation with A^uv ~^ 
0, we can treat the streaking of photoemission in an exact 
fashion. The total wave function |^'(t)) at time t can be 
written as 

|*(r)) = Cd{r)\d)\C)+CUr)\L)\C) 
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FIG. 3: Emission time delays At between satellite and main 
lines for CuBr2, CuCb, and CuF2 with respect to wxuv- In 
the inset, COEs of streaked photoemission spectra within the 
sudden approximation, i.e., V = 0, are given at uxuv = 141 
eV for CuCl2. 
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Dynamics will start by turning on the XUV pulse so 
that the initial state should be |5'(to)) = \ipo) given at 
r = To <C — Txuv- The streaked photoemission spectro- 
gram could be calculated for the main and satellite lines 
as follows: S'„i(k, riR_xuv) = |(k|('0m|*(rmax))P and 
S's(k,TiR_xuv) = KkKV'sl^'lTmax))!^ for a given value 
of Tipt_xuv- Spectrograms are found to converge at 
T = Tmax 3> Txuv- In the actual calculation, we adopt 
To = — 15 a.u. and Tmax = 100 a.u.. 

In the streaked spectrogram 5m(k, tir-xuv) and 
5s(k, TiR_xuv) with distributions of photoelectron ki- 
netic energy, we determine the center of energy (COE). 
COEs of streaked photoemission spectra are presented 
with respect to tir_xuv at wxuv = 141 eV in the right 
panel of Figl^l The metrology of streaked spectra reveals 
a time delay At in the emission of electrons correspond- 
ing to the satellite with respect to those corresponding 
to the main line. This finding is dramatic. Main and 
satellite lines are not simultaneously created, but satel- 
lites are emitted later than main lines by At = 40 as, 73 
as, and 90 as at wxuv = 141 eV for CuF2, CuCl2, and 
CuBr2, respectively. 

The emission time delay At strongly depends on the 
XUV pulse energy wxuv as shown in Figl3] Notably, At 
already disappears at wxuv — > 20 eV, where is 
the binding energy of the satellite, E^ = - Eq k. 131 
eV, which is almost the same for the three compounds, 
as shown in the left panel of Figl^J The same scattering 
potential V also gives rise to extrinsic loss of photoelec- 
tron in photoemission spectroscopy. Extrinsic loss would 
give a change in relative peak strengths or an asymmet- 
ric broadening of peaks observed in the spectroscopy. 
It was argued that extrinsic loss due to extended ex- 
citations in a metal sustains up to the photon energy 
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FIG. 4: (a) Energy-dependent phases of dipole matrix ele- 
ments corresponding to the satellite and main lines of CuBr2, 
CuCl2, and CuF2. (b) Emission time delays Ar between 
satellite and main lines. Photoelectron scattering is taken 
into account within the first-order perturbation theory. 



- 0(1) keV[Ti|, whereas that due to local excitations 
in an electron-correlated insulator disappears at around 
10 - 20 eV above the threshold fl5|. In the inset of FigH 
COEs of streaked photoemission spectra within the sud- 
den approximation, i.e., putting V = 0, at wxuv = 141 
eV for CuCl2 are given, so that At = is in fact con- 
firmed. They indicate that At could be a quantification 
of extrinsic loss contribution to the photoemission spec- 
tra of a given system. 

More figurative understanding of At is possible by not- 
ing that the Wigner-Smith delay is determined by the 
energy derivative of the quantum phase of the matrix el- 
ement. The exact final state I^J^fg-i) in the scattering 
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where E is the energy of the final state, given by i? = 
£k + £-111(8) for the main and satellite, respectively, and 

T — X^k^kCkCk- To the first order in V, we neglect V 
in the denominator. We can then obtain the matrix el- 



ement M„i(s)(k) = (^J^^^g^ jAjV'o)- The quantum phase 
of the matrix element will be (5iii(s)(£k) — arg[Af,„(s)(k)], 
the behavior of which is given in FigUt^a). The time 
delay At between the satellite and main lines is cal- 
culated by At = 9(5s(£k)/9£k - 9i5ni(£k)9£kU+<5-B and 
given with respect to wxuv in FigUJ^b), which could be 
compared with Figl3] Suppression of At with wxuv is 
much faster in the first-order perturbation so that At 
disappears at wxuv — E^ > 5 eV. Perturbation the- 
ory cannot give quantitatively good results compared to 
the exact ones in FigsH] and [3] due to the strong scat- 
tering potential, but could give intuitive understanding 
that the Wigner-Smith delay occurs mostly in satellite, 
especially in the low energy region, i.e., 9(5s(ek)/5ek ^ 

i9(Sin(ek)9ek|£k+<5-E- 

A strong electron correlation between the core level 
and the 3d shell plays an essential role in inducing 
the emission time delay between the main and satel- 
lite lines in copper dihalides. The electron correlation 
U leads to two final states with different screening of 
the core hole created after photoemission. The larger 
the U, the greater the difference. When keeping the 
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FIG. 5: Emission time delays Ar with respect to U at cjxuv ~ 
141 eV obtained by keeping the same initial state by fixing 
the value of Sd + U (= 2.4 eV as in CuCb). All other material 
parameters are taken from CuCb. 

same initial state by fixing the value of Sd + U, as il- 
lustrated in Fig[TJa), we may have |V'm) \d) and 
\4>s) \L) in a limit of large U ^ oo, and |k)|i/)ni) will 
not be scattered by the scattering potential V, i.e., per- 
fect screening, while |k)|'0s) is maximally scattered, i.e., 
no screening. This will cause an increase of Ar, which 
is clearly demonstrated in FiglS] It should also be noted 
that the correlation U scales the potential strength, i.e., 
U = — Vsc(0)[l^. From Ae « U under an assumption 
of small hybridization (See FiglIJa)), main and satellite 
lines would overlap at C/ f« 5.3 eV and Ar would then 
vanish. 



Finally, we note that the transport delay occurring in- 
side the solid was not considered here. However, such 
transport delay could be separated from the Wigner- 
Smith delay because it stems simply from different group 
velocities between satellite and main lines without addi- 
tional scattering [18]. Further, if an atom-thick adsorbate 
is technically possible with copper dihalides, the trans- 
port delay may in principle be removed in the experi- 
ment. 

To summarize, we explored the attosecond streaking of 
3s core lines of copper dihalides by incorporating a model 
of copper dihalide cluster to describe the photoelectron 
scattering. We found that in the emission timing, the 
satellite comes later than the main line. Such an emis- 
sion time delay originates from the electron correlation 
between the core level and the 3d shell, which leads to 
different core-hole screening depending on satellite and 
main lines. The time delay strongly depends on wxuv 
and its disappearance is energetically rapid, i.e., at just 
^ 20 eV above the threshold. To our best knowledge, 
this is the first study to explore the emission timing of 
electrons in a true solid with strong electron correlation. 
Lastly, we claim that the time delay has the same root 
in the extrinsic loss of photoelectrons and it is a quan- 
tification of extrinsic loss. 
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